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SUMMARY 

The  linearized  supersonic  flow  theory  has  "been  applied  to  the 
calculation  of  the  camber  and  twist  required  for  various  lift  dis¬ 
tributions  on  a  triangular  wing  with  leading  edges  behind  the  Mach 


Several  specific  examples  have  been  derived  and  plots  of  the 
cambered  surfaces,  the  corresponding  lift  distributions,  and  values 
of  the  theoretical  drag  coefficients  are  presented.  Thp  examples 
were  chosen  to  avoid  the  singularities  in  lift  distribution  or 
camber  that  appear  in  previously  given  solutions  and  may  be  useful 
for  practical  applications  since  they  avoid  large  pressure  peaks 
and  extreme  angles  of  twist. 


INTRODUCTION 

The  design  of  cambered  wings  requires  a  knowledge  of  the  rela¬ 
tion. between  the  shape  of  the  lifting  surface  and  the  distribution 
of  lift  over  it.  For  airfoils  of  triangular  plan  form  at  supersonic 
speed,  two  examples  of  this  relation  are  known;  namely,  the  lift 
distribution  for  a  flat  wing  at  an  angle  of  attack  (reference  1)  and 
the  camber  required  to  support  a  uniform  distribution  of  lift  over 
the  surface  (reference  2). 

.The  solution  for  the  flat  wing  shows  an  infinite  concentration 
of  lift  along  the  leading  edges,  while  that  for  the  uniform  lift 
distribution  shows  that  an  extremely  large  angle  of  attack  at  the 
center. line  of  the  triangle  would  be  required.  Thus,  both  solutions 
show  singularities  which  the  designer  may  wish  to  avoid.  It  was 
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suggested  to  the  present  author  by  R.  T.  Jones  that  these  undesir¬ 
able  singularities  in  the  lift  and  camber  distribution  could  be 
avoided,  at  a  particular  lift  coefficient,  by  an  appropriate  choice 
of  the  functions  used  to  represent  the  lift  distributions.  In  the 
present  report,  several  examples  of  such  nonsingular  lift  distribu¬ 
tions  and  their  corresponding  camber  shapes  have  been  derived. 

Following  a  brief  review  of  the  method,  a  series  of  functions 
capable  of  representing  a  variety  of  lift  distributions  is  presented. 
It  is  shown  that  the  desired  property  of  the  lift  distribution  and 
camber  shape  can  be  achieved  by  an  appropriate  relation  between  the 
coefficients  of  this  series.  Numerical  examples  showing  lift 
distributions  and  camber  shapes  which  involve  the  first  few  terms  of 
the  general  series  are  then  given. 


METHOD 


Conical  Flow  Theory 


As  in  the  usual  linearized  theory  (references  3  and  4),  the 
boundary  conditions  are  satisfied  in  the  horizontal  plane  rather 
than  in  the  actual  plane  of  the  wing.  The  coordinate  system,  as 
depicted  in  figure  1,  is  chosen  so  that  the  apex  of  the  triangle  is 
at  the  origin.  The  positive  x  axis  extends  back  in  the  plane  of 
symmetry  of  the  triangle,  the  z  direction  is  vertical,  and  the  plane 
of  the  triangle  is  approximately  horizontal.  The  symbols  used  are 
defined  in  Appendix  A. 

Thg  fractional  perturbation  velocities  are  denoted  by  u,  v,  and 
w  (i.e.,  fractions  of  stream  velocity)  in  the  x,  y,  and  z  directions, 
respectively.  The  difference  in  pressure  coefficient  between  upper 
and  lower  surfaces  is  given  by 


=  4u 


and  the  slope  of  the  surface  by 


The  linearized  equation  is  solved  using  the  method  described  in 
reference  4 .  Solutions  are  found  for  a  Mach  number  of  */£,  but  may 
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be  extended  to  other  Mach  numbers  by  the  use  of  the  veil-known  Prandtl 
transformation. 

In  reference  k  it  is  shown  that  a  special  complex  variable 
exists  such  that  if  f(£)  is  any  analytic  function  of 

u(£)  =  real  part  of  f (£)  and  w(£)  =  real  part  of  ^  • 

is  a  conical  solution  of  the  linearized  equation. 

The  complex  variable  £  represents  a  transformation  of  the 
conical  flow  field  which  causes  no  distortion  of  the  airfoil  coordi¬ 
nates  on  the  z=0  plane  inside  the  Mach  cone.,  but  places  the  Mach 
cone  along  the  real  axis  in  the  £  plane  from  £  =  —  co  to  — 1  and 
from  1  to  «v  The  exact  expression  for  £  is 


i  = 


2e 

1+e2 


where 


€  = 


The  use  of  the  variable  £  results  in  a  simplification  of 
boundary  conditions  on  the  plane  z=0,  since  it  can  be  shown  that 
£  approaches  y  +  iz  as  z  approaches  zero, 
x 

To  yield  the  lifting  case,  f(£)  must  be  so  chosen  that  its  real 
part  will  vanish  everywhere  on  the  real  axis  in  the  £  plane  except 
on  the  wing,  which  extends  from  £  =  -m  to  £  =  +m.  In  addition,  the 
real  part  of  f(£)  must  be  antisymmetric  with  respect  to  the  real  axis 
in  the  £  plane.  If  the  wing  lies  entirely  within  the  Mach  cone, 
f ( £)  must  vanish  at  infinity.  This  latter  requirement  arises 
because  in  the  £  plane  the  Mach  cone  is  to  be  found  at  infinity  in 
any  direction  as  well  as  along  the  real  axis  as  previously  described. 

For  wings  inside  the  Mach  cone,  the  integration  for  w  can  be 
taken  along  the  real  axis  from  £  *  =  1  to  £  *  =  £  making  w  zero 
everywhere  on  the  Mach  cone  if  f(£)  is  zero  at  infinity.  The  term 
£*  is  a  dummy  variable  of  integration  replacing  £  in  the  integrand. 

Thus,  to  represent  the  lifting  surface  and  its  pressure  distri¬ 
bution,  f(£)  must  have  the  following  properties: 

1.  f(£)  must  be  an  analytic  function  of  £  . 


x+  */x2-y2-z; 
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2.  f  (£)  must  vanish  at  infinity. 

3.  The  real  part  of  f (£)  along  the  real  axis  must  represent 
the  desired  lift  distribution  of  the  wing. 

4.  The  real  part  of  f(£)  must  he  antisymmetric  with  respect  to 
the  real  axis  in  the  £  plane. 

When  checking  a  prospective  f(£ )  against  these  four  require¬ 
ments,  t  must  he  considered  a  complex  variable.  However,  once 
these  requirements  are  met,  all  subsequent  calculations  are  made  on 
the  z=0  plane,  which  corresponds  to  the  real  axis  in  the  £  plane, 
so  that  t  is  equal  to  y/x,  a  real  variable.  Thus  in  the  final 
expressions  for  u(£),  v(£),  and  £  (£)»  t  ifi  to  y/x* 

An  example  of  a  function  satisfying  the  foregoing  requirement  is 


Here 


tit)  = 


•J m2 — £2 


’  TO 


for  1 £  |  "S  m 


u(^)  =  0  for  |^|  t  m 


u(x,y,0)  = 


oyer  the  yrlrxg. 


Since 


£ 

v(£)  =  real  part  ^-i  J  ■  df»^ 


it  is  found  that 


v(£)  =  -  i  E  ( Vl-my)  for  U  I  =  m 
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The  integration  is  taken  from  x«  »  y|  to  x*  -  x  so  that  the 
Xem  surface  in  which  the  wing  lies  may  at  J-O  on  th e 

S  cone.  Thus,  the  shape  of  the  surface  is  obtained  in  the  form 

z  =  z(x,y). 


Since  the  wing  surface  is  conical,  its  shape  can  he  expressed 
in  the  form  |  =  §  (C),  that  is,  |  =  function  of  £  where 

^  _  Z,  Furthermore,  the  integration  for  z  can  he  accomplished 
in  terms  of  £,  thus 


The  magnitude  of  camber  is  related  to  the  lift  coefficient  by 
finding  the  common  factor  of  u,  w,  and  —  which  will  yield  the 
desired  lift  coefficient. 

total  lift  _  =  1  r  s  =  it  r  r  Uds 

1  q  X  wing  area  S  Jg  j  q  S  Jg  J 


If  u  =  Bu*  where  B  is  the  desired  factor. 


For  the  symmetrical  triangular  wing  extending  from  x=0  to  x=l 
al  nng  the  center  line  and  from  y  =  —  mtoy  =  +  m  along  the  trail¬ 
ing  edge,  the  integration  can  be  accomplished  with  a  single  integra¬ 
tion  in  terms  of  t .  Since  £  is  constant  in  the  element  of  area 
enclosed  by  a  triangle  the  base  of  which  is  dy  along  the  trailing 
edge  and  the  vertex  of  which  is  at  the  vertex  of  the  wing. 


C 


-m 


and,  for  a  symmetrical  distribution, 
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< 


Vd  m 

=  —  p  u*  (t  )<3-t,  or  B 


m 


X 


U*(^)dC 


(2) 


The  coefficient  of  drag  due  to  lift  is  obtained  as  follows: 


The  integration  over  the  wing  can  again  be  reduced  to  a  single 
integration  in  terms  of  £  . 


cd  =  -  5  r 

Jo 


(3) 


No  additional  term  for  leading--edge  suction  is  required  because 
the  pressure  distributions  chosen  have  no  infinite  pressure  peaks 
along  the  leading  edges. 

SERIES  DEVELOPMENT  FOR  AN  ARBITRARY  LIFT  DISTRIBUTION 
Elements  Suitable  for  Series  Representation 
The  analytic  functions 


fo(  c)  =  i  log  ^ 


w 


and 


2n— 1 


7 


WO  - 


n=l,  2,  3>  etc. 


(5) 
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can  be  superimposed  to  yield  an  arbitrary  lift  distribution  extending 
from  -m  to  +m.  The  function  of  f2n-i  is  similar  to  that  given  in 
reference  5 •  The  exponent  is  restricted  to  odd  integers  so  that  the 
lift  distribution  will  be  an  even  function  of  £,  or  symmetric  with 
respect  to  the  center  line  of  the  triangle.  Such  a  solution  would  be 
represented  by  a  series  of  the  form 

11=1 


The  velocity  components  corresponding  to  the  Individual  terms  are 


u2n-i(0  =  real  part 


(6) 


and 


v2nr-i(0  =  real  part 


(7) 


It  can  be  seen  that  this  family  of  functions  has  the  required 
properties  of  being  analytic,  vanishing  at  infinity,  and  having  real 
parts  which  are  antisymmetric  with  respect  to  the  real  axis  in  the 
{;  plane.  Also  the  uprt-i  terms  are  zero  everywhere  on  the  real  axis 
except  at  |  £  |  =  m. 

As  an  example  of  the  requirement  that  the  real  part  of  f (£)  be 
antisymmetric  with  respect  to  the  real  axis,  it  will  be  shown  that 
i(C— nA2-®2)  meets  this  requirement.  For  this  purpose  C  can  be 
taken  equal  to  UlL  since  the  argument  can  be  restricted  to  the 
region  where  z  approaches  zero.  It  is  immediately  apparent  that 
the  real  part  of  i£  is  anti  symmetric  with  respect  to  the  real  axis. 
That  the  real  part  of  i  «^2-m2  this  symmetry  is  not  obvious  but 

it  can  be  seen  as  follows: 
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Let 


£+m  =  rie*<P1,  £-m  =  r^*^2 


where  ri  is  the  length  of  the  line  joining  the  point,  -m,  on  the 
real  axis  to  the  point,  £  ,  in  the  £  plane;  and  cpi  is  the  angle 
this  line  makes  with  the  positive  direction  along  the  real  axis. 
Similarly,  r2  is  the  length  of  the  line  joining  +m  to  £  and 
is  the  angle  this  line  mates  with  the  positive  real  axis. 

Then 


£2-m2  =  (£+m)  (£-m) 


=  rie 


r2eicp2 


3ga) 

=»  i  */rir2  cos  +  i  sin  (~^^j 

and 

real  part  (ivp-^n2)  =  —  */rir2  sin 


From  the  definitions  of  Ti  and  it  follows  that  the  real 

part  of  1  */\  2-m2  is  always  negative  above  the  real  axis  and  always 
positive  helow.  This  is  seen  by  starting  at  a  point  above  the  real 
axis  with  cj>l  and  cp2  less  than  it  and  following  a  path  which  does 
not  pass  through  that  part  of  the  real  axis  between  +m  and  -m. 
Following  a  path  that  does  pass  between  +m  and  -m  is  prohibited 
because  the  real  part  of  f(£)  will  be  discontinuous  there  if  it  is 
antisymmetric.  The  real  part  of  f(£)  is  not  discontinuous  along  the 
real  axis  outside  of  +m  and  -m  because  it  passes  through  zero  there 


i  V  £2-m2  =  i*/rir2  e 


< 
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Orthogonality  Relations 

Thp  upn-i  terms  are  orthogonal  in  the  interval  -m  to  +m  and 
are  zero  at  the  limits  of  the  interval.  Any  reasonable  even  conical 
lift  distribution  can  be  expanded  in  an  infinite  series  of  Uo  and 
Uen-i* 

Let 

£  =  m  sin  0 

for  |  £  |  ^  m  on  the  real  axis ,  then 

j?2nr-i(m  sin  e)  =  i(sin  6  -  i  cos  0)2nr_1 

=  i1_2n+1(cos  9  +  i  sin  0)2It-1 
_  ei(2n-i)e 

and  finally 

fp-n— i  (m  sin  0)  =  (—l)11-1  [cos(2n-l)0  +  i  sin(2nr-l)e  ](8) 

since 

Ugn— i  (m  sin  0)  =  real  part  f  fan— i  (m  sin  0)] 

uotw i  (m  sin  0)  =  (— l)11-1  cos(2n— 1)0  (9) 


To  establish  the  orthogonality  of  these  functions,  it  is  necessary 

r- 

to  evaluate  /  2  u2n_i  u23>-i  d0: 


/ 


77*  7T 

2u2xt-i  u2 j»_i  d0  =  |  Js  .cos(2n-l)0  cos(2i--l)0  d0 


or 


f‘ 


2 

o,  n  ^  r- 


2  iiait—i  u2r— 1 


ru>  u  T  i'| 
if,  n  =  rJ 


(10) 


Equation  (10)  is  the  required  orthogonality  relation 
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Let  h(£)  be  any  even  function  which  is  zero  at  £  =±  m. 
Expanding  h(m  sin  0)  in  the  Fourier  series. 


00 

h(m  sin  0)  =  ^  ban—!  (- l)n  1  cos(2n— 1)0 


and 


n=i 


£ 

b2n- a  =  (— l)n  ^  h(m  sin  0)  cos(2n— 1)0  d0 


(ID 


The  term  aou©  =  ao*  can  be  added  to  h(m  sin  0)  to  give  a  lift 
distribution  which  is  not  zero  at  £  =±m. 

Then 


00 


u'(m  sin  0)  =  ao«  +  ^  h2n— 1(— l)n  1  cos(2n— 1)0 


(12) 


n=i 


or 


00 

u*(£)  =  a0x  +  ^  b2n_i  usn_i(^) 


(13) 


n=i 


Since  the  real  part  of  ^i  log  f°r  |£  |-  m  an(^ 


terms  usn-i  are  'tiie  real  parts  of  i 


0=4 


2  2 
-m 


2n— l 


it  can  "be 


seen  that  the  f  *(£  )  which  corresponds  to  u*(£)  in  equation  (13) 
is 
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which  in  terms  of  £  =  m  sin  0  is 


,  .  sin  0+1 

f’(m  sin  0)  *  i  aQ  log  — — — —  + 

sin  0-1 


£  w*(-i)n  1 

n=i 


[cos(2i>-l)0 


+  i  sin(2n-l)0] 


(15) 


By  the  use  of  the  orthogonality  relation  (equation  (10))  it  can 
he  shown  that  no  net  lift  is  developed  hy  u^n— i  ( 5)  when  n  is 
greater  than  one. 

The  total  lift  is  proportional  to 


J  yuds  =  u(£)d£ 


then 


/m  n  2 

Uen-i(£)d£  =  /  u21i_i(m  sin  0)  m  cos  0  d0 


n-i  n2 
=  (—1)  m  /  cos(2n— l)0  cos  0  d0 


and 


-ml  11211—1 
o 

mjf ,  n=l 

l  <E)a?-{oUx} 


ui  da 


(16) 
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From  this  relation  the  factor  B,  -which  relates  the  magnitude  of 
camber  to  the  lift  coefficient,  can  be  expressed  in  terms  of  aQ 
and  bi,  the  first  two  coefficients  of  the  Fourier  expansion  of 


Therefore,  if 


then 


00 

u(£)  *  B  j^O71  +  ^  ^211—1  ^211— l(5) 


r- 


(5)<U 


aomit+bx®^ 


jt(4a0+bi) 


Condition  for  a  Finite  Value  of  w(5)  at  5=0 

The  expressions  for  Van-i  and  w0  all  have  the  same  type  of 
infinity  at  £=0,  namely,  cosh  1  By  an  appropriate  relation 
between  the  coefficients  of  the  Fourier  series,  the  resultant 
coefficient  of  cosh-1  ^  may  be  made  to  vanish  and  the  infinity 
in  w(  5)  at  t=0  will  thereby  be  avoided. 

According  to  A.  Busemann  in  reference  3>  the  necessary  and 
sufficient  condition  that  there  be  no  infinity  in  w  at  5=0  is 
that  df=0  at  5=0. 

Differentiating  equation  (15.)  , 


(df)e=0  =  i  2a0+  2  (-l)nr"1(2n-l)b2n_1  J  d0  (5=0  when  0=0) 


Therefore  the  necessary  and  sufficient  condition  that  there  he 
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no  infinity  in  w  at  £=0  is 


2a0+ 


00 

^(-l)n_1(2n-l)b2IM.  =  0 
11=1 


(18) 


This  relation  restricts  the  lift  distribution  to  the  extent 
that  the  value  of  u  at  £=0  must  be  less  than  the  average  value, 
provided  that  no  negative  lift  is  allowed  at  any  part  of  the  wing. 
Also  the  value  of  u  at  £=m  depends  on  the  shape  of  the  distri¬ 
bution  of  u  over  the  whole  wing. 


Examples  of  Functions  Contained  in  the  Series  Development 

Evaluation  of  w(£)  and  £  (?)  becomes  increasingly  complicated 
for  higher  values  of  n.  The  following  is  a  list  of  the  expressions 
for  u(  ? ) ,  w(£),  and  §  (?)  of  the  first  three  terms  of  the  series: 
(They  were  derived  by  the  use  of  equations  (la),  (lb),  (^-),  (5),  (6), 
and  (7)*  A  sample  derivation  is  given  in  Appendix  B.) 


Uq(?)  =  * 
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^C)  -  J  (  -  |t|  cos-h)  +  |  [k(  </5^?)  -  E(  -  i  (-) 

»3<£>=^vF? 


(M  =  -  3  /k2"-  4  (1-C2)3/*  +  E  fsin"1  ^32,  <^i>) 
m  m  m  \  Vl-m2  / 


—  —  F  ( sin"1  -^21,  ^  _  4  i  +  1  cosh-1!  i 

m  \  4^  /  m  m  U 

=  15!  cos-H)  (2+£2)  ^?-3  U!cos"H] 


+  ^+g  E  (  */l-m2  )  —  —  K  (  «/l-ffl2)  +  —  (») 
m  mm 

The  preceding  expression  for  ((;)  is  derived  from  an  integra¬ 
tion  of  only  those  terms  of  w  which  are  variable  over  the  range 
Id  <  m  plus  a  constant  term  which  is  needed  along  the  line,  £= 0. 

For  each  &■,  there  will  he  one  term  from  the  integration  of 

cosh-1  which  is  infinite  at  f,=0.  When  the  individual  solutions 
are  superimposed  in  a  series,  the  cosh  1  ^  terms  are  removed  by 

imposing  equation  (l8).  Terms  of  A  (tj)  which  are  to  be  removed  by 
equation  (l8)  have  been  denoted  by  the  symbol  (»)  in  the  table. 

The  elliptic  integrals  F  and  E  are  complete  and  hence  constant 
for  |£|  ^m.  Also  is  taken  to  be  zero  for  is  I*  m,  since  w 

was  to  be  only  the  real  part  of  the  function  resulting  from  an  inte¬ 
gration. 

An  arbitrary  dihedral  angle  can  be  added  to  the  wing  surface 
along  the  wing  center  line  without  changing  the  lift  distribution  as 
long  as  no  part  of  the  wing  is  far  from  the  z=0  plane.  Thus  a  term, 
a  constant  times  £ ,  can  be  added  to  the  equation  of  the  surface 
without  changing  w(£)  or  u(£).  The  value  of  this  constant  determines 
the  dihedral  angle  along  the  wing  center  line.  It  can  be  seen 
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mathematically  that  a  term  £  equals  constant  times  £  will  be 
proportional  to  and  w  (which  equals  will  be  zero. 

Such  a  term  added  to  the  equation  of  the  wing  surface  has  the 
effect  of  shifting  the  stream  sheet  in  which  the  wing  lies  so 
that  it  does  not  pass  through  z=0  on  the  Mach  cone.  However, 
such  a  departure  is  allowable  under  the  condition  that  no  part  of 
the  wing  surface  sheet  lies  far  from  the  z=0  plane.  This 
represents  an  approximation  which  is  valid  to  the  same  order  of 

magnitude  as  the  practice  of  calculating  ^  in  the  z=0  plane 

rather  than  in  the  wing  surface.  For  example ,  Stewart’s  solution 
for  the  flat  triangle  yields  a  V-shaped  wing-surface  cross  section 
which  is  equivalent  to  a  flat  triangle  when  the  dihedral  angle  is 
removed. 


Examples  of  Specific  Lift  Distributions  . 

Four  cases  have  been  worked  out  using  f 0(£) ,  fx(0  j  an<^  fg(C) 
in  various  combinations  which  comply  with  equation  (18) ,  The 

expressions  for  f(£),  u(t;)>  *(£)>  and  £(£)  are  listed  in  each  case. 
Graphs  are  drawn  of  u  and  ^  for  Ct=1.0,  m=0.577  in  figure  2. 

JL 

For  purposes  of  comparison,  the  flat  triangle  and  constant  load 
triangle  are  also  included  in  the  figure. 

The  wing-surface  shapes  of  figure  2  are  drawn  to  a  square  scale. 
(The  unit  of  length  in  the  horizontal  direction  is  equal  to  the 
vertical  unit  of  length. )  The  z=0  plane  is  the  horizontal  line 
above  the  wing.  The  semispan  of  the  wing  is  in  each  case  equal  to 
0.577  of  the  radius  of  the  Mach  circle.  The  leading  point  of  the 
wing  would  be  at  the  origin  of  the  coordinate  system,  a  distance 
ahead  of  the  plane  depicted  in  the  figure  equal  to  the  radius  of  the 
Mach  circle  in  that  plane. 

The  cases  considered  are  most  easily  evaluated  by  the  use  of 

Yi 

the  following  tabulated  expressions  for  u*  v,  and 
u(£)  *  B  |a0u0  + 

n=x 


Vbai—i  U2xri  (£) 
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v(5)  =  B  |^a0Vo  + 

¥0  -  b  [*%  + 

The  term  will  he  zero  for  n>2  in  the  oases  considered 

here. 

The  coefficient  aQ  may  be  taken  equal  to  1  or  0,  If  aQ  is 

equal  to  0,  b  may  be  taken  aB  1,  and  b_  can  be  chosen  with  a 
1  —1  i  3 

value  such  that  the  cosh  —  term  does  not  appear  in  w*  ( t)  or  the 

^  1 

(“)  term  does  not  appear  in  \  (£).  This  choice  is  equivalent  to 

«A» 

satisfying  equation  (l8). 


^  b2nr-i  ■w2n-i 
n=i 

to] 

n=i 


From  equation  (17) , 


Ql 


The  coefficient  of  drag  due  to  lift  was  found  by  a  graphical 
integration  of  equation  (3) . 

A  term  D  |£  |  1b  added  to  j.(£)  to  make  allowance  for  the 

arbitrary  dihedral  angle  at  the  center  line  of  the  triangle.  For 
example,  D  can  be  given  the  value  required  to  place  the  leading 
edges  at  the  same  vertical  height  as  the  center  line.  This  value 
was  used  in  the  cases  considered  here. 
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Case  I 


-i  (*"5 

y(£)  =  2L  |  n/1~?-  ("cosh.-1  I  7^ I  +  cosh.""1  |i=S£  |)  -  | 

'  ='  2s  l  m  \  |  £+m  I  |  £-m  1/  m 

-  £  J^K  (  -/l-ffl.2)  -  E  (  */l-m2)  j- 

§u>  *  £  (£  — 1  |t^I  -  ¥  I^D 

_  £  ^*/l— £2  —  |£  |  cos-h)  ~  “  j^K  (  Vl-m2)  —  E  (  */l-ffl2 )  —  D  |C  |  ^ 


Cl^1 


m  =  0.577 


D  =  1.486 


Cu  =  0.25 


0 

.05 
.10 
.15 
.20 
.22 
.24 
.26 
.28 
.  .30 
•  35 
.40 
.45 
.50 
.55 
.577 


-0.459 

-.429 

-.4o4 

-.383 

-.366 

-.361 

-.356 

-.353 

-.350 

-.347 

-.346 

-.350 

-.362 

-.384 

-.421 

-.459 


0.182 

.188 

.201 

,228 

.271 

.341 

.500 
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Case  II 


td)  = 


log  +  T-s  (S  -*A2-“2)' 


£-m  3m3 

u(£)  -  h  [” +  3?  _ 


—i 

cosh 

l+m  £ 

+  cosh  1 

1 .-a£ 

' 

£+m 

t-m 

-S^ 


_  JL  (1-t2)3/2  +  —  (4  +  m2)  E  ( ,/l-m2)  -  5m2  K  (,/l-ffl2) 
3m3  3m2 


?«) 


aj-s  (i^cosh-1 

Art  ^  m  \  m 


l+m£ 


£+m 


m 


cosh' 


— x 


l-a£ 


) 


3m' 


,/l-C2  -  Ul  cos 


4 

3m1 


(2+^yi-C2-  3  It  |  cos  ^ 


(4+m2)E(  ,/l-m2)  -  5m2  K(y  1-m2)  -  D  |£  | 


Cl-1 


a  =  0.577 


D  =  22.481 


Cp  =  0.26 


1_ 

.05 

.10 

.15 

.20 

.22 

.24 

.26 

.28 

.30 

.35 

.4o 

.45 

.50 

.55 

.577 


Z 

X 

i 

u 

-0.546 

0 

0.196 

-.489 

.1 

.204 

-.442 

.2 

.224 

— .4o6 

.3 

.253 

-.380 

.4 

.286 

-.373 

.5 

.302 

-367 
-363 
—361 
-.360 
-365 
—  382 
-,4io 
—.449 
—504 
-.546 

.577 

.250 

ti|H 
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Case  III 

r (t ) = - 1  r i  *  a  (?-vF^)3 

jt  L  m  3331  j 

u(U  =  2L  4r  (2^+m2)ym2-f 
«  3m 

w(C  )  =  ^{-  ^5  (1-f  )®/2  +  ^3  (2-m2)  E  (/l^)  -m2  K  ( ,/i^2)]} 


MCA  TN  No.  1816 


21 


Case  IV 


f(?)  - 1 1  [108  - :  <  c  +  51  <e 


3/2 


cosh 


— 1 


1+mt,  —1  l-m£ 

*31  +  00Bh  pr 


2^3/2 


3m' 


(l+m2)  E  (V  1-m^ )  -2m2  K  ( ./l-m2 )  j. 


J  7l-m2  t 

6^+m  -1 

_____ 

l+m£ 

_  cosh"”1 

l-m£ 

m  \ 

ir  UUOXl 

^  m 

i;+m 

m 

t-® 

) 


--  (-A?-  Id  cos"1  5)  [(S+^jy^F-Sld  cos-1£] 

in  3m° 

+  -^3  [(!+m2)  E  (a/  l-m2)  -2m2K  (  a/  l-m2 )  ]  -  D  |£|  1 
3hi  J 


O 

l> 

H 

m  =  0.577 

.D  = 

11.984 

CD  = 

z. 

z 

'  X 

i 

X 

't 

u 

0 

-0.517 

0.28 

-0.357 

0 

0.192 

.05 

-.469 

.30 

-.356 

.1 

.198 

.10 

-.430 

•  35 

-.359 

.2 

.216 

.15 

-.398 

.4o 

-.371 

-3 

.245 

.20 

-.376 

.45 

-.394 

.4 

.281 

.22 

-.369 

.50 

-.428 

.5 

.316 

.24 

-.364 

.55 

-.476 

•  577 

•  333 

.2  6 

-.360 

.577 

-.517 
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EVALUATION  OF  RESULTS 

As  can  be  seen  in  figure  2,  the  theoretical  coefficients  of 
drag  due  to  lift  of  the  camber  shapes  considered  here  are  all  larger 
than  that  of  the  flat  triangle.  In  practice,  however,  it  is  possible 
that  case  II  and  case  IV  would  have  lower  drag  than  the  flat  triangle 
at  some  values  of  lift  coefficient  because  of  the  suppression  of 
large  adverse  pressure  gradients. 

Quite  possibly  case  II  and  case  IV  would  be  capable  of  attain¬ 
ing  a  higher  lift  coefficient  than  the  flat  triangle  because  of 
their  favorable  pressure  gradients,  although  the  flat  triangle 
pressure  distribution  would  be  superimposed  on  their  pressures  at 
lift  coefficients  other  than  that  for  which  they  are  designed. 

Case  III  has  the  most  radical  camber,  which  might  cause 
uncalculated  additions  to  the  drag  as  well  as  flow  separation. 

Theoretically  the  dihedral  angle  in  the  surface  shape  is 
arbitrary.  In  practice,  however,  variation  of  this  parameter  might 
result  in  an  optimum  of  characteristics  at  other  angles  then  those 
chosen  in  figure  2. 

While  the  linear  theory  cannot  be  expected  to  apply  for  values 
of  lift  coefficient  greater  than  about  0.25, the  calculations  were 
made  at  a  lift  coefficient  of  one  to  better  illustrate  the  wing- 
surface  shapes.  This  value  of  lift  coefficient,  0.25,  is  also  an 
approximate  limit  beyond  which  the  assumption  that  no  part  of  the 
wing  lies  far  from  the  z  =  0  plane  is  invalid. 

These  solutions  are  applicable  to  swept— back  or  other  wings 
obtained  by  cutting  out  the  rear  of  the  triangle  along  lines  such 
that  no  part  of  the  remaining  wing  lies  within  the  zone  of  influence 

of  the  parts  removed. 

/ 

Ames  Aeronautical  Laboratory, 
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Moffett  Field,  Calif. 


APPENDIX  A 
Symbols 

common  factor  of  u,  w,  and  z  which  relates  the  magnitude  of 
camber  to  the  lift  coefficient 


B 
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°D 

CL 

D 

E 


€ 

t(0 

F 

K 

m 

M 

A£ 

<1 

q 

s 


coefficient  of  drag  due  to  lift 
lift  coefficient 


constant,  the  value  of  which  determines  the  arbitrary 
dihedral  angle  at  the  center  line  of  the  triangle 

elliptic  integral  of  the  second  kind  (either  complete  of 
incomplete) 

y+iz 

complex  variable  equal  to  — ■ — - — . - 

x+  «/x2-y2-z2 


analytic  function  of  the  complex  variable  £  such  that  the 
real  part  of  f(£)  equals  u 

incomplete  elliptic  integral  of  the  first  kind 

complete  elliptic  integral  of  the  first  kind 

horizontal  slope  of  the  leading  edge  of  the  wing 

Mach  number  (equal  to  J2  ) 

the  difference  in  pressure  coefficient  between  upper  and 
lower  wing  surface 

dynamic  pressure  (H 
wing  area 


u,v,w  fractional  perturbation  velocities  as  fractions  of  the  stream 
velocity 

uT  u  divided  by  B 

Y  stream  velocity 


x,y,z  Cartesian  coordinates 


i 


equal  to 


y 

— ,  also  used  as  complex  variable  equal  to 
z 


26 

1+e2 


£*,xf  dummy  variables  used  in  integrations  when  a  limit  of  integra¬ 
tion  is  t  or  x 


2k 
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J  |  absolute  value 

(oo)  terms  in  the  equation  of  the  wing  surface  which  are  to  be 

removed  by  the  superposition  of  other  independent  solutions 


APPENDIX  B 

SAMPLE  CALCULATION  OF  EQUATIONS  IN  TABLE  OF  u,  v,  AND  ^ 


From  equation  (5) 


Since  the  subsequent  calculations  are  made  in  the  plane 
z  -  0,  £  is  a  real  variable  equal  to 


dfi 


-i(i _ L=) 

\m  m*/£  2-m2  ' 


From  equation  (7) 

Vi  «  real  part 


1  _  rs*  l - ZTZ 


^  g  ^  1  P ^  A-~^  2-dM 

mJ1  ^  X'  2-m2 


Vl .  i  p  vid 

BWl  £’ 

P  d£»  =  «/l—  £2—  cosh-1 

/i 


,  _ _ L  dg»  =  r ^  ^ri- 

x  J1  VC2-®2  J  m  */i  ,2~™2 


(Bl) 


(B2) 


/^SELat'-r^-dt  for  t<m 


(B3) 


The  last  equality  holds  because  the  second  term  on  the  right 
of  equation  (B3)  is  imaginary  for  £<m 
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With  the  substitution. 


r2 


1-m2 


_ rgdr _ 

,/i-r2  Vl— ( 1-m2 )  r2 


l-mg)r2  dr 

ji-** 


dr 


V  1— r2  vMl-m2)!-2 


The  two  terms  on  the  right  of  the  last  equality  are  recognized  as 
complete  elliptic  integrals  of  the  second  and  first  kind,  respectively, 
so  that 


A 


d£»  =  EO/L-m2)  -  kC^/l^d2)  for 


(b4) 


The  same  substitution  leads  to  the  result 


These  expressions  all  apply  for  £  negative  as  well  as 
positive  because  w*  is  an  even  function  about  £  =0. 


Thus  from  (Bl),  (B2),  (B3),  (b4),  and  (B5) 


Wi 


-i{^- 


cosh' 


-! 1 


E  CyS^)  +  KCyi-m2)  |  for  |  C  |  <* 
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-H 


Jl— —  cosh  1  —  1  (  sin  -1 


_i  ./I 


J  }  for  |{|J: 

v  yi^  y  J 


From  equation  (lb) 


A 


I  ^ I  vi  ( £ 1 ) 


Since  the  dihedral  angle  at  the  center  of  the  triangle  is 
arbitrary,  the  elliptic  integrals,  which  are  complete  and  hence 
constant  for  It  is  m,  need  not  be  included  in  the  integration. 
However,  to  achieve  the  correct  angle  of  attack  at  the  center  line, 
they  must  be  taken  into  account.  This  operation  can  be  accomplished 
most  simply  by  the  use  of  equation  (la),  which  is 


■VII i  * (£).*■ 


Thus  along  the  center  line,  the  term  of  zi  from  the  elliptic 
integrals  in  the  expression  for  wi  is 


K(  yd?)  -  E(  yi-B?)  dx»  =  J  K-eJ  (] 

The  terms  resulting  from  the  integration  of  cosh  1  |t-[  will 
be  removed  by  a  superposition  of  other  solutions  according  to 
equation  (l8)  to  avoid  an  infinite  angle  of  attack  at  the  center 
line.  Such  terms  are  denoted  by  the  symbol  (®°)  in  the  expression 


*  2^  t  . . . 

The  only  term  of  wi  not  yet  considered  is  —  J  1—  £  .  Inserting 
this  term  in  (B6)  yields 


-id 


i.  V1.  i—  dp  =  i- L/l-£2  -  Id  cos  1  £ 
m  ^2  m  L/ 
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Thus  from  (B6),  (B7)>  and  (B8) 


</hF  -  \i\  cos-1?  +i  k(  )  - 1( yin?  ) 
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